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By using Density Matrix Renormalization Group (DMRG) technique we study the phase diagram 
of ID extended anisotropic Heisenberg model. We analyze the static correlation functions for the 
spin operators both in- and out-of-plane and classify the zero-temperature phases by the range of 
their correlations. On clusters of 64, 100, 200, 300 sites with open boundary conditions we isolate the 
boundary effects and make finite-size scaling of our results. Apart from the ferromagnetic phase, we 
identify two gapless spin-liquid phases and two massive ones. Based on our phase diagram and on 
estimates for the coupling constants known from literature, we classify the ground states of several 
edge-sharing materials. 

PACS numbers: 75.10.Jm, 75.30.Kz, 75.10.Pq, 75.40.Mg 



I. INTRODUCTION 

Recently, an increasing attention has been paid to 
the materials containing edge-sharing Cu02 chains (e.g. 
Rb2Cu2Mo30i2 as in Ref. [l|, NaCu202 as in Ref. or 
LiCuV04 as in Ref. ^ . It has been found by using several 
complementary experimental techniques that the low- 
energy physics in such materials is one-dimensional^. It 
has been also concluded that in such insulating materials 
the spins, localized on the copper ions, interact via ferro- 
magnetic interaction with their nearest neighbors, while 
a considerable next-nearest neighbor (NNN) interaction 
has been argued to be antiferromagnetic. In addition, 
at least in one material (LiCuV04), about 6% exchange 
anisotropy has been measured by using paramagnetic res- 
onance^. Another manifestation of such anisotropy is the 
dependence of the saturation value of the external mag- 
netic field on its direction'^. 

It is widely accepted that to study such systems, a ID 
extended anisotropic Heisenberg model with ferromag- 
netic (F) nearest-neighbor (NN) interaction and antifer- 
romagnetic (AF) NNN one should be used. In order to 
estimate the values of the exchange interactions in these 
materials, temperature and magnetic field dependencies 
of the integrated quantities, such as susceptibility and 
magnetization, have been compared with those of vari- 
ous ID spin models, calculated exactly on small clusters. 
The F-AF Heisenberg model appeared to be the only 
compatible^.. Excluding the symmetry breaking in XY- 
plane, the resulting Hamiltonian could have at most four 
parameters: two in-plane interaction constants and two 
out-of-plane ones. Let us denote these as follows J±, J^, 
Jz and J'z . Such large number of parameters makes it dif- 
ficult to explore the full phase diagram, which would be 
four-dimensional. There exist several parametrizations 
of this system, which use less parameters. For instance, 
one can restrict oneself only to isotropic case (Jz = Jj_, 
J'z = J'j_), ending up with only two parameters, each 
of which may be ferromagnetic or antiferromagnetic (see 
i.e. Ref.Q). Contrarily, one could allow for the same level 



of anisotropy both in nearest and next-nearest neighbor 
channels, which results in three parameters (Jz, J'z and 
J±/ Jz = J±/ J'z) as in Refs.0. Yet, another parametriza- 
tion way, adopted in the present paper, consists of let- 
ting the anisotropy only in the nearest neighbor channel, 
while leaving the NNN interaction isotropic, as in Ref. H. 
This parametrization amounts to have three parameters: 
Jz, J± and J' = J'z = J'j_. Accordingly, the ID extended 
anisotropic Heisenberg model reads as follows (contrarily 
to Ref. S we invert the sign of Jz). 

H = -J.Y.s'^Ui + J±J2'^s:s:^, + sfsf^^ 

i i 

+ j'^S,S,+2. (1.1) 

i 

It is worth noting that the Hamiltonian (|l.ip with 
J± > can be easily mapped onto the one with J± < 0. 
Indeed, for a bipartite lattice, the transformation: 

S-f = (-l)*S'f , where a = x,y. (1.2) 

inverts the sign in front of J± and changes an in-plane 
correlation function (S'f5°^„) by the prefactor (— 1)". 
We choose AF sign of J± to be compatible with our pre- 
vious work^, although in the edge-sharing Cu02 mate- 
rials J± should be negative. Finally, it is worth noting 
that the model (jl.ip can be also represented as a two-leg 
zigzag ladder as show in Fig. [TJ 

From the theoretical point of view, F-AF Heisenberg 
model is a challenge both for analytical and numerical 
methods. The model (|l.ip has been extensively stud- 
ied in its antiferromagnetic region [Jz < 0). For what 
concerns the ferromagnetic region [Jz > 0), a few re- 
sults should be mentioned. First of all, it is worth not- 
ing that despite the non-integrability of the model, there 
exist two points in its parameter space where the an- 
alytic expression for the ground state energy and wave 
functions are known. These are the points: {Jz — — 1, 
J± = Jz, J' = Jz/'^) as proven in Ref. [l^ and {Jz = 
1, J±_ ^ Jz, J' = Jz/'i) as found in Ref. [ill Sec- 
ondly, the famous Haldane conjectur o^^d'^ states that 
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FIG. 1; Mapping of the model p.ip on a two-leg zigzag ladder. 



one-dimensional spin chains, composed of half-integer 
spins, could only have massless excitations and power- 
law-decaying correlation functions. Another conjecture 
- the so-called "ladder conjecture" - comes from nu- 
merical methods^, bosonizatio n^^'^^ and subsequently 
from experiments on a series of ladder materials like 
Sr(n-i)/2Cu(n+i)/20„ (u = 3, 5 . . This conjecture 

states: "spin-i ladders composed of an even number of 
chains have gapped excitations, while those with an odd 
number of chains have gapless excitations" . Neverthe- 
less, the question of whether the model (jl.ip could have 
a gapped ground state and whether the gap can be ob- 
served numerically is still controversial. The existence 
of an astronomically small gap^^ has been predicted by 
means of the Renormalization Group analysis of the ef- 
fective field theory in the limit of two AF chains coupled 
by a weak F inter-chain coupling. A gapped dimer phase, 
surrounded by the gapless spin-liquid ones, has also been 
identified by using the level-crossing analysis of the ex- 
cited states obtained by Lanczos diagonalization on small 
rings22i2i. Similar results were found by using the Quan- 
tum Renormalization Group^^. In addition, a dimerized 
gapped phase^"^ and incommensurate power-law correla- 
tions^i have been found on the AF side of the model 
{Jz < 0). On the other hand, in Ref.'25' by means of 
the perturbation theory around the mean-field solution 
of weakly anisotropic F-AF Heisenberg model, only spin- 
liquid gapless phases have been found. This result might 
be an artefact of the perturbation theory, since a pertur- 
bation can not change the underlying mean-field ground 
state unless one is summing up an infinite number of di- 
agrams. 

Despite the existence of exact Bethe-ansatz solution, it 
is still extremely difficult to obtain closed analytic expres- 
sions for the correlation functions (see e.g^t^) even for 
the isotropic antiferromagnetic Heisenberg model (J' — 
0, J± = Jz). Asymptotic analytical results come mainly 
from the Quantum Field Theory. In Ref. [2^ it has been 
shown, based on the Conformal Field Theory, that the 
long-range behavior (r — > oo) for the isotropic antiferro- 
magnetic Heisenberg model has a highly non-trivial form: 



where a, (3 = x, y, z. As soon as the rotational invari- 
ance is broken, however, the correlation functions of the 
nearest-neighbors anisotropic Heisenberg (XXZ) model 
assume a simple power-law form: 

(S^S^) ^ {-IfA^r'^ + B.,r^+^'^ (1.4) 
(5o^5,^) ^ i-iyAzr^''^ + Bzr-\ (1.5) 

Here {SlSf) = {S^S'f), 77 = - arccos( J,/ J^)/?!^ < 
\Jz/Ji.\ < 1, Ax has been determined in Ref. [2^ and 
Az, Bx, Bz are not known in general for arbitrary values 
of the ratio Jz/ J±- 

In the past we have already considered the ferromag- 
netic {Jz > 0) model (|l.ip in connection with the ques- 
tion of ergodicity of the system's dynamics^ '^°i'^^ , by 
means of the Lanczos and Exact Diagonalization tech- 
niques. Therein a zero-temperature phase diagram has 
been constructed, based on whether the dynamics of the 
^-projection of local spin was ergodic or not. Different 
types of ergodic phases have been identified depending 
on how the non-ergodic constant was approaching zero 
within the finite-size scaling. Another insight into the 
physics of (jl.ip has been made from the entanglement 
studies of the phase transitions in the vicinity of the 
ferromagnetic phase. Two different types of behavior 
can be identified, depending on whether J' or J± is in- 
creased22i^i2i. By means of Lanczos, it was impossible 
to understand deeper the nature of the underlying phases 
without analyzing the asymptotic behavior of the corre- 
lations at large distances and accurate finite-size scaling 
for large clusters. That is why we revisit the phase dia- 
gram of (jl.ip in the present manuscript using the Density 
Matrix Renormalization Group (DMRG). 

This article is organized as follows: after the defini- 
tion of the model and the method used to analyze it in 
Section [ni we briefly revisit the classical phase diagram 
of ([Til) at T = in Section [ml In Section Ml we de- 
scribe the various phases found in the quantum phase 
diagram of (|l.ip . 



II. METHOD OF ANALISYS 

In the present work, we focus our attention on studying 
the spin-i one-dimensional anisotropic Heisenberg model 
with NNN interaction (jl.ip . We use Jz > 0, which cor- 
responds to ferromagnetic coupling. 

We obtain the ground state properties of the Hamilto- 
nian (jl.ip numerically by means of the DMR G^^'^^ tech- 
nique on chains with L = 64, 100, 200, 300 sites, subject 
to open boundary conditions (OBC). In the DMRG cal- 
culations, we maintain up to to = 200 lowest eigenstates 
of the reduced density matrix in the basis of each DMRG 



1.5 - 



- 






- 


q=7l / 






Spiral ± incommensurate 




q=arccos(Jj^/4J') 


Ferro z 




q=0 





0.2 



1.5 - 



Spin-Liquid-I 
(XXZ) 



M-II 



0.5 - Ferro 



'A M-I J$ Spin-Liquid-II 



0.4 0.6 
J7J, 



0.8 



0.2 



0.4 0.6 
J7J, 



0.8 



FIG. 2: Quantum phase diagram of the Hamiltonian (|l.ip (right panel) as compared to the classical one (left panel). The 
dashed regions are the phase boundaries determined with systematic error owing to the finiteness of the cluster size. 



block. The real-space spin-spin correlation functions are 
calculated by means of the finite-system algorithm. 



A few words should be said about the determination 
of the phase boundaries on a finite-size cluster. When 
the correlation length becomes comparable to the cluster 
size, it is not possible to determine the exact location 
and the order of the transition, but only the region in 
the phase diagram, where the transition should occur in 
the bulk. Such region usually shrinks increasing the clus- 
ter size. Moreover, the properties of a cluster with OBC 
converge slower to the thermodynamic limit in compar- 
ison to those of the clusters with e.g. periodic bound- 
ary conditions. On the other hand, the method of level- 
crossing^2^, successfully applied to the determination 
of the phase boundaries in Lanczos and Exact Diagonal- 
ization techniques, cannot be used efficiently in DMRG 
calculations, since clusters with OBC lack many of the 
important symmetries. Nevertheless, DMRG has the im- 
portant advantage of being able to simulate systems as 
large as several hundreds of sites. In an OBC system, 
the boundary effects penetrate inside the system at some 
finite length A;, (Ah/a < 20 in our analysis, where a is 
the lattice constant). Therefore, on a cluster of hundred 
or several hundreds of sites, the central part of the clus- 
ter behaves effectively as a bulk system. The long-range 
part of the correlations can be safely observed in this 
part of the cluster. In the present paper, the transitions 
we deal with are the transitions between a massive and a 
massless phases. In such a case, the phase boundary can 
be estimated from the massive side as the point where 
the correlation length becomes of the order of the cluster 
size, while from the massless side as the point where the 
power-law behavior ceases to be observed. The difference 
between these two points is a measure of uncertainty in 
determining the exact position of the transition line. This 
uncertainty decreases with the system size. 



III. CLASSICAL PHASE DIAGRAM 

Before considering the phase diagram of p.ip it is in- 
structive to take a look at the classical limit of this model 
and its phase diagram at T = 0. The calculations are 
rather straightforward and will be just sketched here. In 
classical case, the spins are represented by the classical 
vectors of constant length s. One has to minimize the 
energy functional of the system, containing the scalar 
products between the nearest and next-nearest vectors, 
subject to the constraint that the length of each vector 
be s. Fourier transform brings the Hamiltonian in the 
"diagonal" form (i.e. the operators S'^ and 5^,, where 
a = x,y, z, are not coupled unless q' = ^l)'- 

H = - J.^cos qS^Stq + cos 2gS,S_, 

q 1 

-l-Ji^cos<z(5^-5^,-H5,^5^,) = (3.1) 

iMsis^^, + Y.^dq){s'^s^-, + sysi^). 

9 q 

Here ^z{q) and 7±(<z) are defined as follows: 

72 (q) = —Jz cos q + J' cos 2q (3.2) 
~ J±cosq + J' cos2q. (3.3) 

In such circumstances it is easy to see that the global 
minimum of the energy is realized by taking only one 
component with q = qo and maximal value Sqg = sy/N 
and putting all the others to zero, go is chosen such that 
either of the following conditions takes place: 

1z{qo) min; Iziqo) < 7±{q),yq 
l±{qo) ^ min; -/^{qo) < jz{q).yq 

In the former case all spins are aligned along the z-axis, 
while in the latter they are parallel to the XY-plane. 
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FIG. 3: (Color online). 300 site cluster, J' — O.2J2, J± — l.GJ^. Example of static correlation functions for in-plane (left 
panel) and out-of-plane (right panel) channels plotted in logarithmic scale in the Spin-liquid-I (XXZ) phase. 
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FIG. 4: (Color online). Correlation functions in the Spin-Liquid-I phase at J' = 0.2 J^, J± = 1.6 Jz for different system sizes in 
the in-plane channel (left panel) and out-of-plane channel (right panel). 




FIG. 5: (Color online). 100 site cluster. Spin-liquid-I phase: critical exponent of the in-plane correlations as a function of J' 
(left panel) and that of J± (right panel). The black curve labelled "exact" is the bosonization result from formula (|1.4p 



5 



The above conditions give rise to the classical phase di- 
agram depicted in Fig. [2] (left panel). In the absence 
of quantum correlations the system possesses a perfect 
long-range order, either commensurate or incommensu- 
rate, with periodicity vectors q as shown in Fig.[2l There 
are three phases in the classical limit of the model (jl.ip in 
the range J' > 0, > 0. The first one is the ferromag- 
netic phase (Ferro) with all spins fully polarized along 
the 2:-axis. The second phase is an in-plane Neel antifer- 
romagnet, while the third one is the in-plane incommen- 
surate spirally ordered phase with the periodicity vector 
q = arccos( J^/4J'). The Ferro phase is separated from 
the Neel one by the line J± = Jz, J' < Jz/^, while the 
Neel phase is separated from the Spiral one by the line 
J_L = 4 J', J± > Jz- Finally, the Spiral phas e is separated 
from the Ferro one by the curve J± = ^8J'(1 — 2J'), 
Jz/i <J'< Jz/2. 



IV. QUANTUM PHASE DIAGRAM 

Quantum fluctuations radically modify the classical 
picture. Long-range correlations with constant val- 
ues of correlation functions, independent on distance, 
are substituted by either quasi-long-range ones with 
power-law behavior at large distances, or by short-range 
exponentially-decaying ones. In the quantum phase dia- 
gram of (jl.ip we observe five phases as shown in the right 
panel of FiglS] 

In our previous studies^ we have located the position of 
one of these phases (Ferro) within the range Q < J± < Jz 
and J' < O.SIJ^. Contrarily to the isotropic Heisenberg 
ferromagnet, where the ground state is L-|-l-times degen- 
erate, in the anisotropic model p.ip the ferromagnetic 
phase has only two degenerate ground states with all 
spins either "up" or "down" . Therefore, the Ferro phase 
is structurally identical to that of the classical phase di- 
agram, the only difference being the shape of the phase 
boundary. 

We now describe the remaining four phases of quantum 
phase diagram in Fig. [2l which we have found by use of 
the methods reported in Section HIl 



A. Spin-liquid-I {XXZ-\\\ie) phase 

Above the ferromagnetic phase ( Jj^ > Jz) and for mod- 
erate values of J' (J' < i(J^/4 — Jz)) we find a phase 
that we called Spin-liquid-I or XXZ-\\]^e because of its 
similarity to the ground state of AF XX Z model. A typ- 
ical behavior of the correlations in this phase is shown in 
Fig. [3l The relevant in-plane correlations are antiferro- 
magnetic and power-law decaying, with periodicity vec- 
tor TT and critical exponent rj^, which is a function of 
both and J': 

(^."^.\<i>-(-l)'rf""- 



For what concerns the out-of-plane channel, the relevant 
correlations are always ferromagnetic: 

with the exponent close to —2. This fact confirms once 
again the analogy with XXZ model. 

Since DMRG deals with open boundary conditions 
(OBC), it is important to keep the boundary effects under 
control. With OBC, a two-site correlation function de- 
pends not only on the distance between the two sites, but 
also on the position of these sites. To minimize bound- 
ary effects, in the measurements of the two-site correla- 
tions, we choose these two sites as symmetric as possible 
with respect to the center of the cluster. In this way, 
the boundary effects, owing to the nonequivalence of the 
lattice points under OBC, can be overcome. We have 
studied the dependence of our results on the cluster size 
for a typical point in the Spin-liquid-I phase. As shown in 
Fig. m a progressive increase of L increases the portion of 
points, lying on a straight line, common for all values of 
L under investigation. The greater is the absolute value 
of r]j_{T]z), the smaller are the boundary effects, since for 
large a exponent the correlations decay fast enough on 
the scale of the cluster size L. 

The dependence of r]± on coupling constants along sev- 
eral lines with fixed values of J±_ and J' are shown in 
Fig. [5l left and right panels respectively. 77 appears to be 
monotonically decreasing as a function of J'. In the case 
when J± is constant, upon approaching the transition 
towards the massive phase M-II (see Fig. 1]± shows a 
tendency to diverge to minus infinity. Peculiarly, all the 
curves cross at the same point around J' « OSJz, which 
means that at this point ri± is independent of J±. 

In order to check the quality of our DMRG data, we 
compare our results for r]± at the line J' = with 
the known "exact" result of bosonization for the XXZ 
model given by the formula (|1.4p . One can see from the 
right panel of Fig[5] that our points fall exactly on the 
bosonization line, meaning that our calculations (DMRG 
and the chosen cluster size) can access the bulk low- 
energy physics, described by the bosonization. Between 
J' = 0.277^ and J' = 0.3Jz there is a change of convex- 
ity of »7±(J_l). At J' fa 0.28Jz, r]±{J±) is approximately 
independent on J±, especially for J± far from 1, confirm- 
ing the presence of a crossing point. The limiting value 
at ^ 00 of the critical exponent for the in-plane cor- 
relations appears to be —1/2, as follows also from (|1.4p . 
This means that for J± large enough J' becomes always 
irrelevant. 



B. Spin-liquid-II phase 

We identify another spin-liquid phase in the range 
J' > OASJz, J± < 1.8J' - 0.2Jz I.e. when J' becomes 
dominant, while the NN interactions (parameterized by 
J_L and Jz) become marginal. In Fig. [5] we report a typ- 
ical correlation picture in the Spin-liquid-II phase, while 
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FIG. 6: (Color online). 300 site cluster, J' — 0.8 Jz Ji_ = 0.5 Jz- Example of static correlation functions for in-plane (left panel) 
and out-of-plane (right panel) channels plotted in logarithmic scale in the Spin-liquid-II phase. 
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FIG. 7: (Color online). Correlation lengths of the intra-chain even-distance correlation functions in Spin-Liquid-II phase at 
J' = 0.5 Jz, Ji_ = 0.8 Jz for different system sizes in the in-plane channel (left panel) and out-of-plane channel (right panel). 
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FIG. 8: (Color online). 100 site cluster. Spin-liquid-II phase: intra-chain critical exponent of the out-of-plane correlations r/z 
(left panel) and that of the in-plane ones rj± (right panel) as a function of J' for several lines with J± =const. 
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the raw data for both in-plane and out-of-plane corre- 
lation functions are shown in the respective insets. It 
is more convenient to plot the correlation functions sep- 
arating the inter- and intra-chain parts. Both in- and 
out-of-plane correlation functions have antifcrromagnetic 
character in the intra- and inter-chain channels with pe- 
riodicity vector q ~ 7r/2. The power-law decay of the 
correlations becomes clear after applying the logarithmic 
scale on both axes. Although the dominant NNN inter- 
action (J') is isotropic, the corrections coming from the 
anisotropic NN ones (J^, and Jj_) induce anisotropy in 
the correlations. 

The out-of-plane correlations appear to have the same 
exponent rj^ for both intra- and inter-chain correla- 
tions. Moreover, for the intra-chain correlations, there 
is an additional modulation by a harmonic term ~ (1 -|- 
asin(27ra;/4)), where x is the distance along a leg of 
the ladder and a <^ 1. Such modulation, might be a 
trace of a less- relevant correlations (i.e. of a correlation 
with smaller exponent) which we are not able to mea- 
sure directly. Contrarily to the out-of-plane channel, in 
the in-plane one the exponents of intra- and inter- 

chain (?y±2) correlations appear to be slightly different, 
as shown in Fig. [HI The cluster-size dependence of the 
intra-chain correlations is shown in Fig. [T] Once again, 
by increasing the system size L, the power-law charac- 
ter of the correlations becomes more and more evident. 
The convergence to the thermodynamic limit, however, 
is somewhat slower in comparison with the Spin-Liquid-I 
phase (see FigjT]) , since now we are measuring the correla- 
tions along a leg of the ladder and the maximum distance 
along a leg is only half of the system size. 

The critical exponents for the intra- and inter-chain 
correlations are not universal and depend on the particu- 
lar values of J' and J± . J' dependence of the intra-chain 
correlation exponents for in-plane rj^^ and out-of-plane 
rjz is shown in Fig. [H] For J± < 0.9, rj^ < —1 and 
r]± > —1, while for J± > 0.9Jz % > —1 and t]± < — 1. 
When J' increases, both exponents tend to —1 as in 
this case the system resembles more and more the XX Z 
model on two noninteracting legs. 



C. Massive phases 

The two massive phases (M-I, M-II) exhibit qualita- 
tively different correlations and the transition line be- 
tween them goes from the point (J' = Jz/4, J± = Jz) to 
the point {J' = OAIJ^, J± = 0.5Jz), as shown in Fig. [51 
The transition lines are determined by the behavior of 
the correlation functions. In a massless phase, the corre- 
lation functions are linear in the log-log scale (power-law 
decay), while for a massive phase they are linear in the 
semi- log scale (exponential decay). 

The first massive phase (M-I) is located in the region 
O.SlJz < J' < 0.43Jz, < J± < Jz- M-I is adjacent to 
the ferromagnetic phase and shows finite magnetization 
in the vicinity of the transition line, which gradually goes 



to zero for J' > 0.35Jz- In the M-I phase, residual fer- 
romagnetic correlations induce a very complex structure 
of the static spin form factor in the z-direction. In addi- 
tion to the residual finite magnetization, which shifts up 
the whole plot of {S^S^_^_^), the boundary effects are en- 
hanced with respect to what we find in the other phases 
and in the in-plane channel. That is why we explore this 
phase on clusters of the maximal reachable size on our 
machines: 300 sites. At short distances {d < 20), in the 
out-of-plane channel the correlations are dominated by 
an exponential incommensurate contribution, as shown 
in the left panel of Fig.O Qualitatively different behavior 
appears for d > 20. At medium distances (20 < d < 250) 
we observe incommensurate long-range correlations with 
apparently no decay. Finally, when the spins are situ- 
ated at the opposite sides of the cluster, the correlations 
between them grow as the spins approach the borders 
of the cluster, being this behavior a clear evidence of a 
boundary effect. 

As it can be seen from Fig. [9l the in-plane correla- 
tions present exponential decay. As we have already men- 
tioned, the model (|l.ip can be also considered as a zig- 
zag ladder. Hence, there could be an a priori distinction 
between intra- and inter-chain correlations. The correla- 
tions we measured distinguish whether the two sites stay 
on the same leg or not. The inter-chain in-plane correla- 
tions have smaller correlation length with respect to the 
intra-chain ones, and go gradually to zero as ^ 0, 
since in that case the two legs of the ladder remain cou- 
pled only in the z-channel. In the left panel of Fig. [TOl we 
plot the J' dependence of the in-plane intra-chain corre- 
lation length for several values of J±. £,±{J') appears 
to be independent of J± at least close to the transition 
line Ferro— >M-I. On a ladder of 100 sites, the maximum 
distance along the legs is of 50 lattice constants. On 
approaching the transition M-I^Spin-liquid-II phase for 
J± = 0.0 and 0.5 J^, reaches its maximum value, while 
for the transition along the line J± = 0.7 Jz, it does not 
occur, since there is a transition towards another massive 
phase (M-II). 

The issue of finite magnetization in the M-I phase, 
rather unexpected, deserves a separate study through the 
finite-size scaling in order to conclude whether it might 
be a bulk feature. In the M-I phase and for the values of 
L accessible to us, we find a close relation between the fi- 
nite residual magnetization and the correlation length 
in the in-plane channel, as shown in the left and middle 
panels of Fig. [TT] Namely, the decrease of the magne- 
tization per site is accompanied by the increase of 
The former appears to be almost size independent for 
J' < OAJz- At J' > OAJz, finite magnetization steps of 
height 1/L appear. This occurs because, for a finite clus- 
ter, the z-component of the total spin can only increment 
by one. The disappearance of magnetization is accompa- 
nied by the divergence of the correlation length in the 
in-plane channel (left and middle panels of Fig. [TTjl . On 
increasing J', reaches a maximal value (of the order 
of L/2) at J' = 0.42 OAAJz and, for greater values of 
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FIG. 9: (Color online). 300 site cluster. A sample of the correlations behavior in the M-I phase. Left panel: distance 
dependence of the out-of-plane static spin form factor. Short-range part of the correlations is shown in the inset. The shift 
due to finite magnetization m — 0.06 is equal to — 0.0036 for the given values of J' and J±. Right panel: semi-log plot of 
the in-plane spin form-factor. Intra- and inter-chain correlations have slightly different correlation lengths. Boundary effects 
manifest starting from approx. d > 250. Lines are guides to the eye. 





0.32 



0.36 0.4 
J7J, 



0.44 0.4 0.5 0.6 0.7 0.8 0.9 1 
J7J, 



implies radically different types of orders in the two mass- 
less phases. 

The in-plane correlation function shows an exponen- 
tial decay with a correlation length = 43.29 for 
J' = 0.6 and Ji_ = 1.2J^). Once multiplied by e''/«^, 
as shown in the right panel of Fig. [T^ the correlation 
function presents incommensurate oscillations as a func- 
tion of d with periodicity vector q {q = 0.14 for J' = 0.6 
and J± = 1.2 Jz) and a phase depending on mod(c?, 4) as 
follows. The whole set of distances d = \i — j\ between 
any two sites i and j on a ID cluster splits into four 
sequences: 



FIG. 10: (Color online). 100 sites cluster. J' dependence of 
the correlation lengths for the in-plane correlations in the M-I 
phase (left panel) and for the out-of-plane ones in the M-II 
phase (right panel) for several values of J±. Shown in this 
panel are the correlations along the legs (do and d2 sequences 
as explained in the text). Lines are guides to the eye. 



doil) - 4/ + 0, 
di{l) = 4^ + 1, 
d2{l) = 4l + 2, 
dsH) ^41 + 3, 



0,4,8,12,. 
1,5,9,13,. 
2,6,10,14, 
3,7,11,15, 



J', a power-law fit appears to be more appropriate. 

In the M-II phase, when both J' and J± become large 
enough, the correlations change qualitatively. In this 
phase, the ratio J±/J' is limited from below and above. 
The out-of-plane correlations exhibit an exponential de- 
cay as functions of the distance between any two sites, 
regardless whether they belong to the same or different 
legs (see left panel of Fig. [T^ . We have also investi- 
gated J' dependence of the out-of-plane intra-chain cor- 
relation length for several values of J± in the M-II 
phase, as reported in the right panel of Fig. [TOl ^z{J') 
grows monotonically from the left border of the phase 
(transition Spin-liquid-I— >M-II) to the right one (transi- 
tion M-II— s-Spin-liquid-II). Such non-symmetric behavior where (pd — ■|mod(c?, 4). 



where Z = 0, . . . , 24 for a 100 site cluster. One can easily 
verify by looking at Fig. [T]that the sequence do{l) cor- 
responds to the case when both i and j belong to the 
same leg and i — j is even. Similarly, the sequence d2{l) 
comprises all the cases when both i and j belong to the 
same leg and i — j is odd. The remaining two sequences 
di{l) and d2,{l) are realized when i and j belong to dif- 
ferent legs. All four sequences appear to be modulated 
by the same harmonic term, but with different phases. 
Peculiarly, the phase shift between do(0 ^^'^ ^2(0 equals 
to TT and the same does the phase shift between di{l) and 
d^{l). Finally, the phases of dQ{l) and di{l) differ by 7r/2. 
Summarizing, we can conclude that the in-plane correla- 
tions of the M-II phase assume the following asymptotic 
form: 



{SlSl^^) - exp(-d/a) cos(gd + 



(4.1) 
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FIG. 11: (Color online), m - magnetization per site (left panel) and - correlation length of the in-plane correlations (middle 
panel) as a function of J' at J±_ = 0.5Jz in the M-I phase for different system sizes. Right panel; - correlation length of the 
out-of-plane correlations in the M-II phase as a function of J' for different system sizes. 
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FIG. 12: (Color online). 300 site cluster. A sample of the correlations behavior in the M-II phase. Left panel: semi-log plot of the 
out-of-plane spin form-factor. Both intra- and inter-chain correlations fit the same line. For d > 100, \{S?S^+a) \ < 10"^ which 
is less than the accuracy of our DMRG method. Right panel: in-plane spin form-factor multiplied by the exponential prefactor 
The points sequences fit quite well on the cosine functions with different phase shifts and with the incommensurate 
periodicity vector q — 0.14. Only distances where the exponential-decaying behavior is established are shown. Lines are guides 
to the eye. 



On a cluster of 100 sites, the out-of-plane correlation 
length ^2 approaches its saturation value of about ^ 
35a, while the system undergoes the phase transition to- 
wards the Spin-Liquid-II phase, and hence £.z/L ^ 0.35 
for L — 100. This value should be compared to L/2 = 50 
since we are speaking about the intra-chain correlations. 
In the right panel of Fig[TT]we report the J'-dependence 
of for different system sizes. The ratio £,'^/L ranges 
from 0.3 for L = 300 to 0.375 for L = 64 while the "un- 
certainty" range, i.e. the width of the transition on a 
finite system, decreases as expected. 



V. CONCLUSIONS 

We have performed an extensive numerical study of 
the model (jl.ip within the following range of the Hamil- 
tonian parameters: Jz = I, < J± < 2.5, < J' < 1.0. 
We have constructed a zero-temperature phase diagram 
based on the DMRG calculations for 64, 100, 200, 300 site 
system with open boundary conditions. 

We have identified several phases within this phase di- 
agram, based on the behavior at large distances of the 
correlation functions. We have found two phases with 
quasi-long-range behavior, among which one appears to 
be a generalization of the ground state of XX Z model, 
and another possesses rather complex correlation picture, 
which distinguishes between intra- and inter-chain dis- 
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tances. We have discovered two massive phases. One of 
them (M-II) separates the two massless phases, while the 
other (M-I) divides a massless phase from the ferromag- 
netic one. Although these two massive phases present 
qualitatively different correlations, they share a common 
feature. Both of them develop when one of the NN cou- 
plings is less efficient respect to the other and to J' . In- 
deed, in the phase M-I, J± is limited respect to the other 
couplings, while in M-II the relative weight of goes to 
zero with respect to J_l and J'. 

A further insight into the nature of these phases can 
be obtained by studying more specific properties of their 
ground states, like entanglement, or searching for exotic 
order parameters. For what concerns real materials, for 
Rb2Cu2Mo30i2 we have^ J' « 0.37 and from Fig.[l]we 
can conclude that its ground state should be in a massive 
phase (M-I or M-II depending on the anisotropy). The 



same conclusion holds also for Li2ZrCu04 (J' « Q.2&Jz)^ 
Pb2[CuS04(OH)2] (J' « 0.43 J^) and CS2CU2M03O12 
(J' « 0.387^) as estimated in the Ref. [H NaCu202 
was found^ to have J' ~ 2Jz and assuming that the 
anisotropy is not too strong we conclude that its ground 
state is in the Spin-Liquid-II phase. Finally, for LiCuV04 
the estimates for J' range from 2AJz to 3.5 Jz. These val- 
ues, under the assumption of weak anisotropy, bring us 
again to the Spin-liquid-II phase. 
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